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Preface 

The three measurement systems were developed and built by J. Louthain. M. Hayduk 
and R. Bussjager built and characterized the fiber lasers. J. Haus developed the mode-locked 
laser model used to simulate the performance of the fiber lasers. 


1. Introduction 


Mode-locked erbium-doped fiber lasers are inexpensive compact sources of ultrashort 
pulses in the 1.55 pm wavelength regime. Ultrashort pulse sources have many potential 
applications including high-repetition-rate sources for telecommunications, local area networks, 
or RF digitization. Recent advances have been made in photonic analog to digital conversion. [1] 
Many passive mode-locking techniques have been presented in the past operating both in normal 
and anomalous dispersion regimes. Recently semiconductor saturable absorbers in the form of 
bulk, multiple quantum well (MQW) or single quantum well - saturable Bragg reflector (SBR) 
have been used to mode-lock linear fiber cavities [2-7]. Complex Ginzburg-Landau equations 
have also been used to model mode-locked laser operation successfully [3, 5]. In particular, 
Kutz et. al. and Haus et. al. made a direct comparison between simulations and experimental 
results on a fiber laser. [3, 5] The dispersion and nonlinearity of the components can be used to 
accurately model the cavity. 

In this paper we present measurements of components of a fiber laser mode-locked by 
MQW saturable absorber. Our cavity is high loss and its non-linear dynamics are dominated by 
the saturable absorber. The construction of our compact, polarization insensitive, passively 
erbium-doped fiber laser employing a MQW saturable absorber and a chirped fiber Bragg 
grating is shown in Figure 1. This is a high loss cavity with an estimated 60 % loss per round 
trip. The measured typical autocorrelation and optical spectrum of the fiber laser is shown in 
Figure 2. [5,8] 
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Erbium Doped Fiber 



Fig. 1. Schematic diagram of the laser cavity. Outputs " 1" and "2" are used to characterize the laser, but 
the normal output is on the left. 



Fig. 2. Autocorrelation spectrum (3.4 ps pulse width) and optical spectrum (0.92 nm 
spectral width) of output pulses from the laser cavity. 

We measured the dispersion of the fiber and Bragg gratings as well as the non-linearity of 
the MQW saturable absorber. Dispersion of optical fiber has become an essential performance 
factor, as bit rates and spectral bandwidths increase in telecommunication applications. Precise 
dispersion information is also important in designing ultra-fast optics. The nonlinearity of the 
saturable absorber is substantial since the structure contains 50-75 periods of 100 A InGaAs 100 
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A InAlAs layers. This large saturation of absorption nonlinearity was required to overcome the 
high loss of the cavity. Unfortunately, large refractive nonlinearity accompanies the absorption 
nonlinearity. The dispersion in the chirped fiber Bragg grating was used in the laser to 
compensate for the large refractive nonlinearity of the saturable absorber. [8] Precise dispersion 
management is required for stable mode-locked operation. [ 5-7] 

The work performed in the measurement of mode-locked fiber laser components can be 

summarized by three efforts: 

1. We developed a measurement system to characterize the dispersion of optical 
fiber. We use spectral interference to measure the dispersion of fiber lengths of 
less than 1 meter to within 0.1-picosecond-resolution. 

2. We built a measurement system to characterize the dispersion of optical fiber 
chirped Bragg gratings. We measured the dispersion of chirped Bragg gratings 
using a phase-shift measurement technique. 

3. We investigated the nonlinearity of the multiple quantum well saturable 
absorber. The z-scan technique was performed on the saturable absorber to 
determine both the refractive and absorptive nonlinearity. 

2. Dispersion measurements of erbium-doped optical fiber 

We used a balanced Michelson interferometer and a broad-band LED, 60nm full-width- 
half-maximum (FWHM), centered at 1540 nm to measure the dispersion of the actual fiber used 
in the laser cavity. Often times dispersion in a small piece of fiber (< 1 m) can differ 
significantly from the average value of the entire spool. In addition, high-loss erbium and 
ytterbium-doped fiber cannot be measured by conventional time-of-flight measurements. Using 
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this method, we were able to measure the linear dispersion of less than a half meter of fiber down 
to as low as 4 fs/nm. The experimental setup is shown in Figure 3. One arm of the 
interferometer includes the fiber length used within the cavity and the other arm is a variable 
free-space section. The light travels through one port of the 3-dB coupler and reflects off the 
broad-band mirrors in each arm. [9] 

The following interference spectrum was collected at the other port and measured by the 
optical spectrum analyzer 

/ = V + V + 2V,,.cos(«^(^a)), (1) 

where 

(l>{(o) = ico)d - {(o)L , (2) 

f A > 4,- 5 Pjs ■> are the respective intensities and propagation constants for the free space 

and dispersive element arms, co is the optical angular frequency, d is the length of the free 
space arai, and L is the length of the dispersive element arm. [4,9,11] 



Figure 3. Schematic of Michelson interferometer used to measure dispersion of optical components. 
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If we expand in a Taylor series about cd = o)^ corresponding to yl-^o collect all constant 
terms into an unknown phase term , ^(a>) becomes 


^(<») = ^(,+ p^{ol)^)L —(<y®o) +g/^3(^o)^(^ ^o) 

where is a constant, c is the speed of light, p^ is the total group velocity dispersion, and p, is 
the third order dispersion of the fiber under test.[4,9,ll] We neglect terms higher than third 
order. Adjusting the free space arm until p^io}^)L = dlc makes the second term drop out, 
leaving the second order and higher terms. 

The wavelength point of symmetry where the phase shift through the free-space path is 
the same as the fiber path occurs at A = Adjusting the free-space pathlength will tune the 
interferometer to different wavelength symmetry points. For instance, the wide center peak in 
Figure 3 corresponds to this symmetry point where A = /!„. Also, from this adjustment one can 
determine the sign of the dispersion. If increasing the free-space path causes the symmetry point 
to move to higher wavelengths, the p^io),) term will be negative. If the symmetry point moves 


to lower wavelengths, the P 2 (®o) dispersion term will be positive. 

The interference spectrum at the output port contains nulls where = {2n + \)7t and 

peaks where = 2 n 7 t as shown in Figure 4. We can now write Eq. (3) in the following form 
for the peaks in the interference pattern: 


h 


+Lp = 0 

A„ 6 Afi Ag 


where A^ is the symmetry point in wavelength, N is the peak location, and A/^ is the wavelength 
of the corresponding peak.[10] We set up a system of simultaneous equations for the three 
unknowns <p^, p^, and p^. Each equation corresponds to a maximum (or minimum) of the 
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interference pattern spectrum with the corresponding value of N. The dispersion can be 
expressed as delay time dispersion, D , in units of ps/nm-km using the following relationship 


£) = - 



( 5 ) 


where is the group velocity dispersion, c is the speed of light, and is the wavelength. [6] 



Figure 4. Typical interference spectrum measured by the optical spectrum analyzer. The associated phase 
is plotted on the right axis by locating the peaks and valleys and assigning the appropriate phase. 

Using the setup in Figure 3, we measured the dispersion of standard single mode fiber 
and erbium-doped fiber. We used data acquisition and instrument control and a PC running 
appropriate software to solve Eq. (4) automatically. See Appendix I for the Labview^^ program. 
This automation allowed us to take several readings. The dispersion of standard single mode 
fiber is published as approximately 
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ps/nm-km, for 1200 nm <A< 1600 nm. 


( 6 ) 



K 

A' 


where A is the operating wavelength, A^ is the zero dispersion wavelength, and is the zero 
dispersion slope.[13] Typical values are 1301.5nm</lo ^1321.5nm and <0.092 ps/nm-km. 
We used the typical value of 0.90 ps/nm-km for -So.[13] 

Figure 5 shows how our experimental results compare with the published dispersion of 
STANDARD FIBER for a 42.5-cm fiber length. At 1550 nm the dispersion is about 16.9 + 1.4 
ps/nm-km and about 17.6 ± 1.4 ps/nm-km at 1560 nm. Notice how our measurements diverge 
from the published results as we moved away from 1550 nm. Our source was centered at 1540 
nm with a spectral width of 60 nm FWHM. Therefore, these variations are most likely due to 
low power levels in the tails of the LED spectrum.[9] 

Figure 6 shows our experimental results for a 38-cm erbium-doped fiber. The power 
levels in the erbium-doped fiber were lower than in the standard fiber measurements due to the 
absorption that occurred at 1530 nm. This contributed to more fluctuation in the measurements 
which in turn required more data points. Again, these measurements will be more accurate 
around 1550 nm. At 1550 nm the dispersion is about 8.5 ± 1.6 ps/nm-km and at 1560 nm it is 
about 10.4 ± 1.6 ps/nm-km.[9] 
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Figure 5. The squares are individual measurements of the dispersion of standard single mode fiber using 
the spectral interference method. The solid line is a linear fit to the data. The top and bottom dashed 
lines are the published dispersion specifications of standard fiber for a zero dispersion wavelength of 
1301.5 nm and 1321.5 nm, respectively. 



Figure 6. The squares are individual measurements of the dispersion of the erbium-doped fiber using the 
spectral interference method. The solid line is a linear fit to the data. The erbium-doped fiber had a 
concentration of about 1200 parts per million and an outside fiber diameter of about 113 gm. 
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3. Dispersion measurements of chirped fiber Bragg gratings 


The chirped fiber Bragg gratings are much more dispersive (over two orders of 
magnitude) than standard optical fiber. Therefore, a time-of flight type measurement is more 
appropriate than an interferometric technique. Our experimental setup is shown in Figure 4. In 
this technique, laser light is intensity modulated with an RF signal and is sent through an equal 
power splitting 2x2 fiber coupler (3dB coupler) to the grating mirror under test. The light 
reflects off the grating and travels back through the splitter to a high-speed detector. Using an 
HP Network Analyzer, we measure the phase shift between the signal reflected off the grating 
and the input RF signal. As the wavelength of the tunable laser changes the phase shift changes 
due to the dispersion. The phase shift is then converted to group delay versus wavelength to 
illustrate the dispersion. [12] 

The chromatic dispersion, D(A,), is obtained from the difference in the group delay at two 
closely spaced wavelengths. The dispersion is the average delay over the wavelength range. 
The resolution of the dispersion measurement is determined by the spectral width of the source. 
The spectral width of our tunable laser was 0.1 nm full-width-half-maximum (FWHM). 
Therefore, we were not able to measure the fine structure of the grating. The end fibers were 
angle cleaved to reduce spurious reflections from the end facets of the fibers. The following 
relation was used to calculate the chi'omatic dispersion coefficient in ps/nm: 


Z)(2) = 10'' 


[<D(/l,)-<D(2,)] 
[360 X / X A/l] 


( 7 ) 


where /l = (A| + 2j)/2, ,) are the measured phase shifts in degrees, -/l^ is the 

difference in nm between the two optical wavelengths, and / is the RF modulation frequency in 
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l lz.[12J We used a modulation frequency of 100 MHz. The following simplified relation will 


find the absolute group delay in picoseconds 10'“ x [O(/l)]/[360x /].[12] 


RF in 


Network Analyzer 


RF out 



Detector 



Modulator 




Tunable Laser 

Device 
under test 


Q 


3 dB Coupler 

Figure 7. Setup for chromatic dispersion measurement of a reflection fiber Bragg grating. 


We used the differential phase shift technique to measure the dispersion and reflection 
spectrums for three very different gratings. We employed data acquisition and instrument 
control and a PC running appropriate software to calculate the group delay. See Appendix II for 
the Labview'^^ program. Figures 8-10 show the group delays. The dispersion in ps/nm 
corresponds to the slope of the curve. The grating in Figure 8 is an ordinary highly reflective 
Bragg grating with a bandwidth of about 1.5 nm. It is very narrow with a relatively flat 
dispersion over the spectral width. Figure 9 shows the dispersion of a linearly chirped grating 
with a bandwidth of about 5.2 nm. The average slope from 1553 to 1556.5 nm is -10.3 ± 0.2 
ps/nm. Finally, in Figure 10 we show the dispersion of a 10 nm wide Bragg grating. This 
grating was built to be a wide-band reflector without any regard for dispersion characteristics. 
The grating has positive dispersion below 1556 nm and negative dispersion above 1556 nm. Due 
to etalon effects in the tunable laser, each of these plots have oscillations with a period of about 
0.35 nm.[9] 
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Figure 8. Reflection magnitude and group delay of a narrow-band high-reflection Bragg grating. 
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Figure 9. The reflection magnitude and group delay of a wide-band chirped grating (3M CS-98 
3312 serial # 8217-7008). The slope over 1553 to 1556.5 nm is -10.3 ± 0.2 ps/nm. 









W avelength (nm) 

Figure 10. The reflection magnitude and group delay of a wide-band chirped grating (3M CS- 
96-2006 serial # 6229-1006). 

4. Nonlinearity measurements of multiple quantum well saturable 
absorbers 

Next, we measured the nonlinearity of the saturable absorbers using the z-scan technique. 
The nonlinear absorption is the mechanism that provides the mode-locking of the laser. The 
refractive nonlinearity balances the dispersion within the cavity. We were able to determine the 
nonlinear absorption due to saturation of absorption and the refractive nonlinearity using the 
same experimental set-up. We determined the nonlinear absorption from an open aperture z- 
scan. This technique was first described by Sheik-Bahae.[15] 

Our Z-scan set-up is shown in figure 11. In order to get the resolution required we 
focused the laser beam source down to a spot-size of 6 pm. We used a PC running appropriate 
software to automate the data acquisition. See Appendix III for the Labview^^ program. 

The normalized open-aperture transmittance z-scan is 
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( 8 ) 


nz,s 


=i)=I 

A/ -0 


(ot + !)■’'■ 


for |c/o| < 1, where 



a is the linear absorption coefficient, p is the two photon absorption (2PA) nonlinear saturation 
of absorption, is the open aperture peak irradiance, Wj, is the beam waist, L is the sample 
thickness, and X is the wavelength. [15] Using these relationships we determined the nonlinear 
absorption, by numerical evaluation. 



Aperture 
Sample ■ 



-ZM —► z 


Figure 11. Z-scan set-up for measuring saturable absorbers. 

The refractive nonlinear index in saturable absorber material is more difficult to extract 
since the z-scan output is dominated by the saturation of absorption. A thorough numerical 
evaluation done by Sheik-Bahae indicated that the nonlinear refractive index was determined 
within less than 10% uncertainty provided that ^o(O) < 1 and \/3/2k')\ < l.[15] The second 
condition is met when Re( 2 ^^^) is larger than or in other words when the nonlinear 

refraction dominates over the nonlinear saturation of absorption. For cubic nonlinearity the 
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index of refraction is expressed in terms of nonlinear indices H 2 (esu) or y (nr/W) through 
n = n^, +'-~\E\'= n„ +/I, (10) 

where uo is the linear index of refraction, E is the peak electric field (cgs), and / represents the 
irradiance (MKS) of the laser beam within the sample.[15] The refractive nonlinearity can be 
determined by dividing the closed aperture normalized z-scan by the open aperture normalized z- 
scan. The result is a new z-scan where ATp.y agrees within 10% of that obtained from a purely 
refractive z-scan, given that the two conditions mentioned above are satisfied. From this result 
the refractive nonlinearity can be calculated by the following relation 

- ( 11 ) 

0 . 406 ^ 4 ^/ 0 ( 1 - 5 )“'' 

where ATp.y is the difference between the peak and valley closed/open normalized z-scan, S is the 
aperture linear transmittance, k=27i/A, and /o is the on-axis irradiance at the focus.[15] 

In Figure 12, we show a plot of the normalized Z-scan. From the open aperture case, we 
were able to determine the nonlinear absorption, p, of -1.1 cm/W. The solid line is the closed 
aperture case divided by the open aperture case. From this normalized plot, we were able to 
determine the nonlinear refractive index, y, as -4.0x10'* cm^/W. The relation /3/2ky was -3.4 
and therefore did not satisfy the relation required to determine the refractive nonlinearity within 
less than 10% uncertainty. The terms «, and / are related through the following conversion 
formula = {cn^ /40;T)y in mV(watt-sec) or esu, where c is the speed of light in a vacuum. [15] 

In Figure 13, we compared the absolute z-scan transmittance of a number of MQW 
saturable absorbers. The quantum well region of the samples consisted of 50 periods 100 A 
InGaAs wells and 100 A InAlAs barriers on an InP substrate, except #1948 which consisted of 
75 periods of 80 A InGaAs and 100 A InAlAs. The wavelength of the source was 1549 nm. The 
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absolute transmittance is calculated by dividing the optical power transmitted through the sample 
(Ps) divided by the optical power transmitted through free space (Pfs = 2.1 mW). We found the 
transmittance to range from 0.25 to 0.5. The higher the peak the more saturation of absorption. 
Figure 14 shows that sample #1948 had the highest saturation of absorption and sample #1442 
had the lowest. From the plotted normalized transmittance we calculated the absorption 
coefficient, P, of each of the samples. Using equations 8 and 9, we determined Ps of -2.5 cm/W, 
-2.9 cm/W, -1.8 cm/W, -1.1 cm/W, and -1.0 cm/W for sample #s 1948, 1590, 1643, unknovm, 
and 1442, respectively. Figure 15 shows the absolute z-scan transmittance of a 500 pm InP 
substrate. The InP substrate thickness was much larger than the Raleigh range of the laser (zo) 
and therefore one cannot use the thin sample relationships to calculate the nonlinear effects. 
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Figure 12. Z-scan of saturable absorber # 1590. The dashed line is the normalized transmittance for the 
open aperture case and the dotted line is for the closed aperture case (at 50% transmittance or S=0.5). 
From the open aperture case, we were able to determine the nonlinear absorption, P, of-1.7 cm/W. The 
solid line is the closed aperture case divided by the open aperture case. From this normalized plot, we 
were able to determine the nonlinear refractive index, y, as -4.0x10*^ cmVw. 
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Figure 13. This plot shows the absolute transmittance open aperture z-scan of a number of different 
multiple quantum well saturable absorbers. The samples are listed from highest to lowest transmittance 



Figure 14. This plot shows the normalized transmittance of a number of different multiple quantum well 
saturable absorbers 1.00 to 1.35 pm thick. The samples are listed from highest to lowest normalized 
transmittance 
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Figure 15. This is the absolute z-scan transmittance of a 500 )im InP substrate. The InP substrate 
thickness was much larger than the Raleigh range of the laser (zq) and therefore one cannot use the same 
relationships to calculate the nonlinear effects of the substrate alone. 

5. Conclusions 

A mode-locked fiber laser design incorporating a MQW saturable absorber and a chirped 
fiber Bragg grating is reliably mode-locked with pulse widths as short as 2 ps. No polarization 
control in the cavity is required. The dispersion in the cavity and the nonlinearity of the MQW 
saturable absorber was measured. The dispersion and nonlinearity measurements can be used to 
simulate the mode-locked operation of the laser. Our fiber dispersion measurements were 
consistent with published values, but there was some variability due to the stability of the 
interference measurements and the limited spectral range of our source. Koch et. al used a 
super-continuum fiber laser with a much higher signal to noise ratio over a very wide spectrum. 
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allowing for measurements over the entire spectral range from 1.1 pm to 1.7|j,m.[14] Koch et. 
al. also used two broad-band 3dB fused couplers instead of the bulk beam-splitters and mirrors to 
simplify alignment and to increase stability. [14] The measured dispersion of different fiber 
Bragg gratings was consistent with the type of gratings, but there were some fluctuations 
inherent in the tunable laser. These periodic fluctuations in the reflectivity measurements and 
the dispersion measurements were most likely due to the etalon effects within the tunable laser 
source. The 0.1-nm FWHM tunable laser was adequate for our purposes, since the spectral 
bandwidth of the mode-locked fiber laser was about 1 nm FWHM. If finer dispersion 
measurements are required, a more stable and narrower linewidth tunable laser source must be 
used.[9,12] 

The nonlinearity results were consistent with MQW structures determined with different 
methods in other papers.[16] The absorptive nonlinearity can be determined fairly easily using 
this technique. The refractive nonlinearity was much more difficult to determine due to the 
strong absorptive nonlinearity, but the sign of the refractive nonlinearity and the order of 
magnitude was easily determined. Some sources of error in our measurements would be due to 
our lack of an accurate measurement of the spot size of the laser. We estimated the spot size of 
the laser to be approximately 6 pm referenced to the e^ power points. A more accurate 
characterization of the beam and spot size would lend more confidence to the measurements. In 
addition, a measurement of known reference material would allow us to calibrate the z-scan set¬ 
up. There was also some variability in the saturation of absorption depending upon where in the 
x-y plane of the MQW the laser was focused. Finally, there could have been some errors 
because we are measuring a compound structure of the multiple quantum well stack and the InP 
substrate. 
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Appendix 


I Dispersion Calculationjinear.vi 

Front Panel: 


Address | 


Xo Beginning Wavelength 
(microns) 2 

Delta X Between points 
(microns) 


stepper I 
^I^BAddr'essj 

iis. 


Steps/micron 
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This program calculates the dispersion. D in ps/nm/km. of a short length of optical 
fiber at a specific wavelength. It displays the optical interference spectrum centered 
at the wavelength of interest. Then the user selects the peaks and valleys using the 
cursors. The center wavelenth is determined as the midpoint between cursors 
"Val R -1" and "Val L-1." The user pushes the "calculate dispersion" button and the 
program calculates two dispersion values one from the location of the peaks and one 
from the location of the valleys. 
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Dispersion Calculation lincar.vi Program Description: 

This program calculates the dispersion, D in ps/nm/km, of a short length of optical liber at a 
specific wavelength. It displays the optical interference spectrum centered at the wavelength of 
interest. Then the user selects the peaks and valleys using the cursors. The user pushes the 
"calculate dispersion" button and the program calculates two dispersion values one from the 
location of the peaks and one from the location of the valleys. 

Block Diagram: 
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II Sweep Phase Wavelength.vi 

Front Panel 


^ 51550 00 ^ ■^S^RnrTn ^ i —^from tunetole lasaH 


^lo IW " 


frdmtun ablg (aaarj > . 

%i\ ■ start Sweep 
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JgpM i (3PlBadclte$SQf'tutl«ibl^Ja5ar 


error out' ‘ ■ 
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Sweep Phase Wavelength.vi Program Description: 

This VI calculates the group delay of a chirped bragg grating versus wavelength. It controls the 
tunable laser and takes the group delay data from the network analyzer. The network analyzer 
drives the electro-optic modulator and measures the phase difference after the tunable laser light 
has reflected back from the chirped bragg grating. It plots the group delay and sends the data to 
an ASCII file. From the slope of the group delay one can detennine the dispersion or the amount 
of chirp in the grating in (ps/nm) 


Sweep Phase Wavelength.vi Block Diagram: 
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Ill z-scan.vi 
Front Panel: 


1 








Z-scan.vi Program Description: 

This program controls an adjustable stage and takes detector measurements at user specified 
increments. You can set the detectors to read directly or taken a mean/median of N readings. 
The information is sent to an ASCII file for further analysis 
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